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1. Introduction 



Nucleus-nucleus collisions at relativistic energies have been intensely studied over the last 
two decades. The main goal of these efforts is to understand the properties of strongly interacting 
matter under extreme conditions of high energy and baryon densities for which the creation of a 
quark-gluon plasma (QGP) is expected ^ Q]. Fluctuations of physical observables in heavy ion 
collisions (see e.g., the reviews §]) may provide important signals regarding the formation of 
a QGP. Measuring the fluctuations one might observe anomalies of the onset of deconfinement 
[||, ^] and dynamical instabilities when the expanding system goes through the 1-st order transition 
line between the quark-gluon plasma and the hadron gas Furthermore, the QCD critical point 
may be signaled by a characteristic pattern in fluctuations Km. With the large number of particles 
produced in heavy ion collisions at CERN SPS and BNL RHIC energies it has now become feasible 
to study fluctuations on an event-by-event basis. 

Let's consider some region in phase space as an acceptance region for particles produced in 
the collision and measure the multiplicity distribution within this acceptance. Than to quantify 
multiplicity fluctuations one uses the scaled variance: 

where (...) denoted event-by-event averaging. By selecting particles with certain quantum num- 
bers one can study fluctuations of electric charge, strangeness, charm, etc. 

Ratio of two species (e.g K/n) in the acceptance also can be studied. To characterize ratio 
fluctuations one usually uses different measures such as a, V, F, <3>, etc. One can also analyze 
correlations between different species in terms of the correlation coefficient: 

_ (AN A AN B ) 
Pab = -fj2 • U-A) 

((AN A ) 2 ) ((Ayv*) 2 ) 1 

Two different acceptance regions can be taken to study long range correlations. For example 
forward-backward correlation measurements which have been performed by the STAR Collabora- 
tion and will be discussed here. We note that higher moments of a distribution also can be used to 
explore the early stage of the colliding system (see e.g. [^]). 

A powerful tool for studying event-by-event physics is the transport model. It allows to imple- 
ment precisely an experimental acceptance, study centrality and colliding energy dependences on 
final results, etc. Results of two transport models will be shown here: the Hadron-String-Dynamics 
(HSD) [g |(J and Ultra-Relativistic-Quantum-Molecular-Dynamics (UrQMD) [JTTJ, |2|. These 



models provide a rather reliable description (see, e.g. Refs. [13, 14, [Tl]]) for the inclusive spectra 
of charged hadrons in A+A collisions from SIS to RHIC energies. Due to the fact that there is no 
explicit quark-gluon phase in HSD and UrQMD, these models can be used to examine the exper- 
imental data owith respect to the presence of 'new physics'. We mentioned that an explicit phase 
transition from hadronic to partonic degrees of freedom is implemented in Parton-Hadron-String- 



Dynamics (PHSD) [15], which will opens new possibilities in studying heavy-ion collisions. 
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2. Multiplicity Fluctuations 



In each A+A collision only a fraction of all 2A nucleons interact. These are called participant 
nucleons and are denoted as N p roj and N p arg for the projectile and target nuclei, respectively. The 



nucleons, which do not interact, are called the projectile and target spectators, N\ 
mdN' s arg =A 



pro) 
S 



N p p mJ 



■N P arg . 



HSD: 158A GeV | 



CJ3 
■- 




Figure 1: The scaled variance co p rg for the fluctuations of the number of target participants, N p rg . The HSD 
simulations of a> P " g as a function of N p "' J are shown for different colliding nuclei, In+In, S+S, Ne+Ne, O+O 
and C+C at £^=158 AGeV. 



In each sample with N P = const the number of target participants fluctuates around its mean 
value, (N p arg ), with the scaled variance (o' p rg . From an output of the HSD minimum bias simula- 
tions of A+A collisions at 158 AGeV we form the samples of events with fixed values of N p ro} . 
Fig. [I] presents the HSD scaled variances co p ' rg as functions of N p roj for different colliding systems. 
The fluctuations of target participants are quite strong for semi-peripheral collisions. These large 
fluctuations in the number of participants strongly influence the fluctuations of extensive observ- 
ables such as multiplicity [fl7]], electric charge and baryonic number [18]. Experimental data of the 
NA49 Collaboration for the multiplicity fluctuations also show an enhancement for semi-peripheral 
collisions (Fig. |||). Though both transport models HSD and UrQMD show a similar enhancement 
for the full acceptance, this dependence become flat when apply experimental acceptance (which 
is in the forward hemisphere). This effect can be explained by the nucleus-nucleus dynamics. 

The consequences of the asymmetry between projectile and target hemispheres - introduced 
by fixing the number of projectile participants N p ro} - depend on the A+A dynamics. According to 



Ref. [|1_9|] different models of hadron production in relativistic A+A collisions can be divided into 
three limiting groups: transparency (T-), mixing (M-), and reflection (R-) models. The rapidity 
distributions resulting from the T-, M-, and R-models are sketched in Fig. [| 

The HSD model (as well as UrQMD) shows only a small mixing on initial baryon flow and is 
closer to the T-model (cf. [[H]]). This supports the findings from Ref. [13] about the influence of 
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Figure 2: The scaled variance of the multiplicity distribution for negatively (upper panel), positively (mid- 
dle panel) and all (bottom panel) charged particles as a function of the number of projectile participants 
Np R0J compared with model simulations in the NA49 acceptance (HSD and UrQMD predictions were taken 
from jl6t]). The statistical errors are smaller than the symbols (except for the most peripheral points). The 
figure is adapted from jl7||. 



the partonic degrees of freedom on the initial phase dynamics which should increase the mixing by 
additional strong parton-parton interactions. 

To decrease fluctuations in the participant number (which result in observable fluctuations), 
one needs to consider the most central collisions. As seen from Fig. [j] for the most central collisions 
these fluctuations vanish. Moreover one sees good agreement between data and transport models 
for more central collisions in Fig. ^. 

An ambitious experimental program for the search of the QCD critical point has been started 
by the NA61 Collaboration at the SPS [2C]. The program includes collecting very central collisions 
varying in the atomic mass number A and in the collision energy. This allows to scan the phase 
diagram in the plane of temperature T and baryon chemical potential near the critical point as 
argued in Ref. [20]. One expects to 'locate' the position of the critical point by studying its 'fluctu- 
ation signals'. High statistics multiplicity fluctuation data will be taken for p+p, C+C, S+S, In+In, 
and Pb+Pb collisions at bombarding energies of Ei a b=l0, 20, 30, 40, 80, and 158 AGeV. We have 
considered these collision systems within the HSD (as well as UrQMD) transport model [21]. The 
study thus is in full correspondence to the experimental program of the NA61 Collaboration [p0|]. 
A monotonic energy dependence for the multiplicity fluctuations is obtained in the HSD transport 
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Figure 3: The sketch of the rapidity distributions of the baryon number or the particle production sources 
(horizontal rectangles) in nucleus-nucleus collisions resulting from the transparency, mixing and reflection 
models. The spectator nucleons are indicated by the vertical rectangles. In the collisions with a fixed number 
of projectile spectators only matter related to the target shows significant fluctuations (vertical arrows). The 
figure is adapted from [ |l9| | . 
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Figure 4: The HSD results for a) in A + A and p + p collisions for the rapidity y > 0, which is in corre- 
spondence with the NA61 program at SPS [ 20 1 . The 1% most central collisions are selected by choosing the 
largest values of Np WJ (see [^TJ for details). 



model (Fig. ^|). Thus, the expected enhanced fluctuations - attributed to the critical point and phase 
transition - can be observed experimentally on top of a monotonic and smooth 'hadronic back- 
ground' . Our findings should be helpful for the optimal choice of collision systems and collision 
energies for the experimental search of the QCD critical point. 
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3. Ratio Fluctuations 

One of the possible measures to characterize fluctuations in the ratio Rab = Na/Nb is a J7|,|22|]: 

, _ ((ARab) 2 ) 



(Rab) 2 ' 

which can be expanded in the following way (see [ p3| ] for details): 
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Figure 5: The scaled variance a> % of pions (left) and cy y „ (3J) of K/% fluctuations (right) at different cm. 
energies ^/snn in the GCE, CE, and MCE ensembles (dashed lines) as well as from the HSD ttansport model 
(solid lines). 

The experimental data for Na/Nb fluctuations are usually presented in terms of the so called 
dynamical fluctuations 



yn 



sign (a 



1 1/2 



(3-3) 



where a m j x corresponds to the mixed events procedure which is calculated in the same way as 



a ( |3.1| ) but for uncorrected particles from different events. Details about the mixed events proce- 



dure can be found e.g. in ||23J]. 

Fig. HI presents the scaled variance (O n of pions and Od yn (3_3) of K/n fluctuations. Please note 
a different behavior of HSD model results for ft) and a as a function of energy. The scaled variance 
ft) increases with energy while the dynamical fluctuations a decreases. The reason is seen from 
Eq. \2: c 2 0)/(N) where (N) is average multiplicity increasing with energy. The difference 
between transport and statistical models also changes when we go from ft) to a. The difference 



becomes larger with energy for the scaled variance [31] but for a it is more pronounced for smaller 
energies. 

In Fig. |6| the HSD results of Od yn for the K/n, p/n and K/p ratios are shown in comparison 



with the experimental data by the NA49 Collaboration at the SPS [25] and the preliminary data of 
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Figure 6: The HSD results for the excitation function in <7 ( / v „ for the Kj n,K j p, p / '% within the experimental 
acceptance (solid line) in comparison to the experimental data measured by the NA49 Collaboration at 
SPS |Q and by the STAR Collaboration at RHIC @ H ||l . The UrQMD calculations are shown by 
dotted lines. See [||, ^ for details. 



the STAR Collaboration at RHIC [g§ g^, g|, ggj. The available results of UrQMD calculations 
(from Refs. [25, 32, |^]) are also shown by the dashed lines. 

The HSD results presented in Fig. [6| correspond to the centrality selection as in the experiment: 
the NA49 data correspond to the 3.5% most central collisions selected via the veto calorimeter, 
whereas in the STAR experiment the 5% most central events with the highest multiplicities in the 
pseudorapidity range \rj\ < 0.5 have been selected. 
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One sees that the UrQMD model gives practically a constant crf^, which is by about 40% 
smaller than the results from HSD at the lowest SPS energy. This difference between the two 
transport models may be attributed to different realizations of the string and resonance dynamics 
in HSD and UrQMD: in UrQMD the strings decay first to heavy baryonic and mesonic resonances 
which only later on decay to 'light' hadrons such as kaons and pions. In HSD the strings dominantly 
decay directly to 'light' hadrons or the vector mesons p, ft) and K*. Such a 'non-equilibrated' string 
dynamics may lead to stronger fluctuations of the K/% ratio. 

At the SPS energies the HSD simulations lead to negative values of a c i yn for the proton to 
pion ratio. This is in agreement with the NA49 data in Pb+Pb collisions. On the other hand HSD 
gives large positive values of ffj^j at RHIC energies which overestimate the preliminary STAR 
data for Au+Au collisions [28]. For o%f n only preliminary STAR data in Au+Au collisions are 
available [^] which demonstrate a qualitative agreement with the HSD results (Fig. ||). The HSD 
results for ov^ show a weak energy dependence in both SPS and RHIC energy regions. 

An interesting feature is a strong 'jump' between the SPS and RHIC values, seen in the middle 
panel of Fig. |6[ in the HSD calculations which is caused by the different acceptances in the SPS and 
RHIC measurements. The influence of the experimental acceptance is clearly seen at 160 A GeV 
where a switch from the NA49 to the STAR acceptance leads to the jump in af f n by 3% - middle 
panel of Fig. ^. On the other hand, our calculations for Pb+Pb (3.5% central) and for Au+Au (5% 
central) collisions - performed within the NA49 acceptance for both cases at 160 A GeV - shows 
only a very week sensitivity of a v £ on the actual choice of the collision system and centrality - cf. 
the coincident open circle and triangle at 160 A GeV in the middle panel of Fig. |6[ 



4. Forward-Backward Correlations 

Correlations of particles between different regions of rapidity have for a long time been con- 
sidered to be a signature of new physics. The observation of such correlations in A+A collisions 



at RHIC energies by the STAR Collaboration [34, 35] has therefore elicited a lot of theoretical 
interest. 



Fig. [7] shows the dependence of the forward-backward correlation coefficient p (1.2) as a func- 



tion of ri gap on the bin size and centrality definition within the HSD model (see [36] for details). 
The dependence of p on r\ gap is almost flat, reflecting a boost-invariant distribution of particles 
created by string breaking in HSD. The right top panel of Fig. [7| demonstrates also a comparison of 
the HSD results with the STAR data [|4|, 35]. One observes that the HSD results exceed systemat- 
ically the STAR data. However, the main qualitative features of the STAR data - an approximate 
independence of the width of the pseudo-rapidity gap r\ gap and a strong increase of p with cen- 
trality - are fully reproduced by the HSD simulations. Note that choosing smaller centrality bins 
leads to weaker forward-backward correlations, a less pronounced centrality dependence, and a 
stronger dependence on the bin definition. The physical origin for this is demonstrated in Fig. ||. 
As the bin size becomes comparable to the width of the correlation band between Np and N'j/ , 
the systematic deviations of different centrality selections become dominant: the same centrality 
bins defined by Np and by N"ff contain different events and may give rather different values for the 
forward-backward correlation coefficient p . 
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Figure 7: The HSD results for the forward-backward correlation coefficient p for 10% (top) and 2% (bottom) 
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Figure 8: The histogram shows the distribution of HSD events with fixed number of participating nucleons 
Np and fixed reference charge particle multiplicity N'jf . The same centrality class (20-22% as an example) 
- defined in various ways - contains different events. 



Thus, the experimental analysis for different bin sizes and centrality definitions - as performed 
here - may serve as a diagnostic tool for an origin of the observed correlations. A strong specific 
dependence of the correlations on bin size and centrality definition would signify their geometrical 
origin! 
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5. Conclusions 

Our analysis has shown that the fluctuations in the number of participants strongly influence 
observed multiplicity fluctuations. To avoid them one should consider the most central collisions 
with rigid events selection. 

We have considered C+C, S+S, In+In, and Pb+Pb nuclear collisions from Ei a b= 10, 20, 30, 40, 
80, 158 AGeV, which is in full correspondence with the experimental program SHINE started at 
SPS. A monotonic energy dependence for the multiplicity fluctuations is obtained in the HSD trans- 
port model. Thus, the expected enhanced fluctuations - attributed to the critical point and phase 
transition - can be observed experimentally on top of a monotonic and smooth 'hadronic back- 
ground'. Our findings should be helpful for the optimal choice of collision systems and collision 
energies for the experimental search of the QCD critical point. 

It has been found that the HSD model can qualitatively reproduce the measured excitation 
function for the K/n ratio fluctuations in central A+A collisions from low SPS up to top RHIC 
energies. Accounting for the experimental acceptance as well as the centrality selection has a 
relatively small influence on Od yn and does not change the shape of the Od yn excitation function. 
The HSD results for a% also appear to be close to the NA49 data at the SPS. On the other hand a 
comparison of the HSD results with preliminary STAR data in Au+Au collisions at RHIC energies 



is not fully conclusive: Od yn from HSD calculations is approximately in agreement with data [29] 
for the kaon to proton ratio, but overestimates the experimental results [^|] for the proton to pion 
ratio. New data on event-by-event fluctuations in Au+Au at RHIC energies will help to clarify the 
situation. 

The forward-backward correlations have been studied within the HSD transport model for 
= 200 GeV. It has been shown that strong forward-backward correlations arise due to an 
averaging over many different events that belong to one 10% centrality bin. In contrast to average 
multiplicities, the resulting fluctuations and correlations depend strongly on the specific centrality 
trigger. When the size of the bins decreases, the contribution of 'geometrical' fluctuations should 
lead to weaker forward-backward correlations and to a less pronounced centrality dependence. 
Note, that the 'geometrical' fluctuations discussed here are in fact present in all dynamical models 
of nucleus-nucleus collisions. Thus, they should be carefully accounted for before any discussion 
of new physical phenomena is addressed. 
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